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We discuss the decay of axion walls bounded by strings and present numerical simulations of the decay process. 
In these simulations, the decay happens immediately, in a time scale of order the light travel time, and the average 
energy of the radiated axions is {ua) — 7ma for Valrua — 500. (uJa) is found to increase approximately linearly 
with ln{va/ma)- Extrapolation of this behavionr yields (uia) — GOrria in axion models of interest. 


1. Introduction 


The axion is the quasi-Nambu-Goldstone boson 
associated with the spontaneous breaking of the 
Upq{1) symmetry which Peccei and Quinn pos¬ 
tulated 1^. Its zero temperature mass is given 
by: 


TOa ~ 6 • 10“®eV • N ■ 


IQi^GeV 


( 1 ) 


where Va is the magnitude of the vacuum expec¬ 
tation value that breaks Upq{1) and is a pos¬ 
itive integer which describes the color anomaly 
of Upq{1). The axion owes its mass to non- 
perturbative QCD effects. At temperatures high 
compared to QCD scale, these effects are sup¬ 
pressed and the axion mass is negligible. There is 
a critical time ti, such that ma{ti)ti = 1, when 
the axion mass effectively turns on. The corre¬ 
sponding temperature Ti ~ 1 GeV. 

In case that inflation occurs with reheat tem¬ 
perature higher than Va, axion strings form and 
radiate axions till ti. At U each string becomes 
the boundary of N domain walls. If A^ = 1, the 
network of walls bounded by strings is unstable 
and decays away. There are in this case three con¬ 
tributions to the axion cosmological energy den¬ 
sity: axions that were radiated by axion strings 
before ti, axions from vacuum realignment, and 
axions that were produced in the decay of walls 
bounded by strings after U. Here, we discuss the 
third of these contributions. 

That axions walls bounded by string decay 
predominantly into barely relativistic axions was 


suggested in Ref. §• D. Lyth 1^ also discussed 
the wall decay contribution and emphasized the 
uncertainties affecting it. Nagasawa and Ka- 
sawaki [Q performed computer simulations of the 
decay of walls bounded by string and obtained 
{uja)/ma — 3 for the average energy of the radi¬ 
ated axions. Our simulations, presented in sec¬ 
tion 2, are similar to those of Ref. |Q] but done 
on larger lattices and for a wider variety of ini¬ 
tial conditions. We obtain {iuja)lma — 7 where 
Vairfia — 500. We attribute the difference be¬ 
tween our result and that of Ref. [|| to the fact 
that we give angular momentum to the collapsing 
walls. 

The lattice sizes which are amenable to present 
day computers are at any rate small compared to 
what one would ideally wish. The lattice constant 
must be smaller than b = (-s/Xua)”^ to resolve 
the string core. On the other hand the lattice 
size must be larger than rn~^ to contain at least 
one wall. Hence the lattice size in units of the 
lattice constant must be of order ^^4 x or 

Triad Triad 

larger if the simulations are done in 2 dimensions. 
Present day computers allow lattice sizes of order 
4000 X 4000, i.e. ~ 100. In axion mod¬ 
els of interest ^ = 10^®. To rem- 

rria 10 ^eV 

edy this shortcoming, we study the dependence of 
{uja)/ma upon the and find that it increases 

approximately as the logarithm of that quantity; 
see section 2. 
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The following toy model incorporates the qual¬ 
itative features of interest: 

+r]Va(t)i 

= + ^d^(j)2d^4i2 

(2) 

where (/) = </>! -|- is a complex scalar field. 
When rj = 0, this model has a C/(l) global symme¬ 
try under which (j){x) e^°‘(j){x). It is analogous 

to the Upq{1) symmetry of Peccei and Quinn. It 
is spontaneously broken by the vacuum expecta¬ 
tion value {(j)) = Uae*“. The associated Nambu- 
Goldstone boson is the axion. The last term in 
Eq. @ represents the non-perturbative QCD ef¬ 
fects that give the axion its mass rua- We have 
rria = ./fj for Va » Wa- 

When rj = 0, the model has global string solu¬ 
tions. A straight global string along the z-axis is 
the static configuration: 

(/)(x) = Ua/(p)e*® (3) 


where (z, p, 9) are cylindrical coordinates and 
f{p) is a function which goes from zero at p = 0 
to one at p = oo over a distance scale of order 
S = The energy per unit length of the 

global string is 

/i ~ 27ry pdp^\^4>\'^ = Trvlln{VXvaL) (4) 

where L is a long-distance cutoff. For a network 
of strings with random directions, as would be 
present in the early universe, L is of order the 
distance between strings. 

When T] 0, the model has domain walls. If 
Va ^ rUa, the low energy effective Lagrangian is : 


= -duad^a + mlvl 


cos-1 

Va 


( 5 ) 


Its equation of motion has the well-known Sine- 
Gordon soliton solution : 


a{y) 

Va 


2-Kn -\- 4 tan ^ exp{niay) 


( 6 ) 


where y is the coordinate perpendicular to the 
wall and n is any integer. Eq. (|^ describes a do¬ 
main wall which interpolates between neighboring 
vacua in the low energy effective theory (^. 



Figure 1. A small cross-section of the universe near 
the QCD phase transition, in A = 1 axion models. 
Lines are domain walls, open (filled) circles are up¬ 
ward (downward) going strings. 

The thickness of the wall is of order m~^. The 
wall energy per unit surface is cr = SruaVa in the 
toy model. In reality the structure of axion do¬ 
main walls is more complicated than in the toy 
model, mainly because not only the axion field 
but also the neutral pion field is spatially varying 
inside the wall |^. When this is taken into ac¬ 
count, the (zero temperature) wall energy/surface 
is found to be: 

a ~ 4.2 f^m^fa ^ 9 ruafl (V 

with fa = Va/N. 

In the early universe, the strings are stuck in 
the plasma and are stretched by the Hubble ex¬ 
pansion. However with time the plasma becomes 
more dilute and below temperature [D 

the strings move freely and decay efficiently into 
axions. Because this decay mechanism is very ef¬ 
ficient, there is approximately one string per hori- 
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zon from temperature T* to temperature Ti ~ 1 
GeV when the axion acquires mass. 

A cross-section of a finite but statistically sig¬ 
nificant volume of the universe near time ti was 
simulated with the assumption that the axion 
field is randomly oriented from one horizon to the 
next just before the axion mass turns on at time 
ti 1^,^. Fig. 1 shows the result of this simulation. 
The average density of walls at time ti predicted 
is approximately 2/3 per horizon volume at time 
ti- 


2. Computer Simulations 

We have carried out an extensive program 
of 2D numerical simulations of domain walls 
bounded by strings. The Lagrangian (2) in finite 
difference form is 


L = 




i{n,t)] + 


+ (</2(n + J,t) - </2(n,t 

4 ^ 


+ {(j)2in,t)f - 1 


+ - 1 ) 


(9) 



Figure 2. Speed of the string core as a function of 
time for the case l/rUa = 400, y/XIma = 10, and 
w = 0. 



Figure 3. Position of the string or anti-string core as 
a function of time for l/nia = 1000, A = 0.0002, D = 
2896, and v = 0. The string and anti-string cores 
have opposite z. They go through each other and 
oscillate with decreasing amplitude. 

where n labels the sites In these units Ua = 1, 
the wall thickness is l/nia = l/i/i? and the core 
size is (5 = 1/VX. The lattice constant is the 
unit of length. In the continuum limit, the dy¬ 
namics depends upon a single critical parame¬ 
ter, niaS = malVX- Large two-dimensional grids 
4000 X 4000) were initialized with a straight 
domain wall initially at rest or with angular mo¬ 
mentum. The initial domain wall was obtained 
by overrelaxation starting with the Sine-Gordon 
ansatz (j)i + i(j )2 = exp(itan“^exp(ma?/)) inside 
a strip of length D between the string and anti¬ 
string and the true vacuum (t/i = 1 , 1 / 2 = 0) out¬ 
side. The string and antistring core were approx¬ 
imated by (pi + i(p 2 = —tanh(.583 r/(5) exp(=F*6*) 
where (r, 0) are polar coordinates about the core 
center, and held fixed during relaxation. Stable 
domain walls were obtained for l/(^ma) ^3. A 
first-order in time and second order in space al¬ 
gorithm was used for evolution with a time step 
dt = 0.2. The boundary conditions were periodic 
and the total energy was conserved to better than 
1%. If the angular momentum was nonzero, then 
the time derivative p was obtained by a finite dif¬ 
ference over a small time step. 

The evolution of the domain wall was stud¬ 
ied for various values of -s/A/ma, the initial wall 
length D and the initial velocity v of the strings in 










4 



Figure 4. Decay of a wall at successive time intervals 
At = 1.2/ma for the case nia^ = 100, A = 0.01, 
D = 524, and v = 0.6. 



'i, 'i, 

(d (f) 


Figure 5. Energy spectrum at successive time in¬ 
tervals At = l/rria, with Wmax = Vs -I- mi and 
Iw = 181n(ui/ma)/ln(u)max/ma) -|-2, for the case 
= 500, A = 0.0032, D = 2096 and v = 0.25. 


the direction transverse to the wall. Fig. || shows 
the longitudinal velocity of the core as a function 
of time for the case m~^ = 400,-s/A/ma = 10 
and u = 0. An important feature is the Lorentz 
contraction of the core. For reduced core sizes 
<5/7 s ^ 5, where js is the Lorentz factor associ¬ 
ated with the speed of the string core, there is 
’scraping’ of the core on the lattice accompanied 
by emission of spurious high frequency radiation. 
This artificial friction eventually balances the wall 
tension and leads to a terminal velocity. In our 
simulations we always ensured being in the con¬ 
tinuum limit. 

For a domain wall without rotation (?; = 0), the 
string cores strike head-on and go through each 
other. Several oscillations of decreasing mag¬ 
nitude generally occur before annihilation. For 
7 s ~ 1 the string and anti-string coalesce and 
annihilate one another. For 7 s ^ 2, the strings 
go through each other and regenerate another 


wall of reduced length. The relative oscillation 
amplitude decreases with decreasing collision ve¬ 
locity. Fig. H shows the core position as a func¬ 
tion of time for the case m~^ = 1000, 'JXjma = 
14.3, £1 = 2896. 

We also investigated the more generic case of 
a domain wall with angular momentum. The 
strings are similarly accelerated by the wall but 
string and anti-string cores miss each other. The 
field is displayed in Fig. ^ at six time steps for the 
case m~^ = 100, A = 0.01, D = 524 and v = 0.6. 
No oscillation is observed. All energy is converted 
into axion radiation during a single collapse. 

We performed spectrum analysis of the en¬ 
ergy stored in the (j) field as a function of time 
using standard Fourier techniques. The two- 
dimensional Fourier transform is defined by 




2i7r 


Px'^x 


Py1 


fin)m 
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where V = L^Ly. The dispersion law is: 


L fr, “^TTPy 

i 2 I 2 — cos —-cos ——^ 


+ m2. (11) 


Fig. I shows the power spectrum dE/d In u) of the 
(j) field at various times during the decay of a ro¬ 
tating domain wall for the case m~^ = 500, A = 
0.0032, D — 2096, and v = 0.25. Initially, the 
spectrum is dominated by small wavevectors, k ^ 
•nia- Such a spectrum is characteristic of the do¬ 
main wall. As the wall accelerates the string, the 
spectrum hardens until it becomes roughly 1/k 
with a long wavelength cutoff of order the wall 
thickness l/rua and a short wavelength cutoff of 
order the reduced core size 6/js- Such a spectrum 
is characteristic of the moving string. 



Figure 6. (w) as a function of time for l/rUa = 500, 
D = 2096, V = 0.25 and A = 0.0004 (solid), 0.0008 
(long dash), 0.0016 (short dash), 0.0032 (dot) and 
0.0064 (dot dash). After the wall has decayed, (lj) is 
the average energy of radiated axions. 

Fig. H shows the time evolution of {ui)jma for 
m~^ = 500, u = 0.25,13 = 2096 and various val¬ 
ues of A. By definition, 

( 12 ) 

_ _ 

P P 

where Ep is the gradient and kinetic energy stored 
in mode p of the field (j). After the domain wall 
has decayed into axions, (w) = (wo). Fig. shows 
the time evolution of {uj)/ma for m~^ = 500, A = 
0.0016, D = 2096 and various values of v. 



Figure 7. Same as in Fig. except A = 0.0016, and 
V = 0.15 (solid), 0.25 (long dash), 0.4 (short dash) 
and 0.6 (dot). 

When the angular momentum is low and the 
core size is big, the strings have one or more os¬ 
cillations. In that case, {uja)/ma — 4, which is 
consistent with Ref. jij. We believe this regime 
to be less relevant to wall decay in the early uni¬ 
verse because it seems unlikely that the angular 
momentum of a wall at the QCD epoch could be 
small enough for the strings to oscillate. 

In the more generic case when no oscillations 
occur, {oja)lrna — 7. Moreover, we find that 
{uJa)lma depends upon the critical parameter 
\f\lmai increasing approximately as the loga¬ 
rithm of that quantity; see Fig. 6. This is con¬ 
sistent with the time evolution of the energy 
spectrum, described in Fig. 4. For a domain 
wall, (w) ~ mo. whereas for a moving string 
(w) ~ ma ln(-\/Ana 7 s/mQ). Since we find the de¬ 
cay to proceed in two steps: 1) the wall energy 
is converted into string kinetic energy, and 2) the 
strings annihilate without qualitative change in 
the spectrum, the average energy of radiated ax¬ 
ions is (Wq) ~ TOa ln(-\/Aua/ ruo). Assuming this 
is a correct description of the decay process for 
'/XVa! mo 10^®, then {oja)lma 60 in axion 
models of interest. 

3. Conclusions. 

We have presented results of our computer sim¬ 
ulations of the motion and decay of walls bounded 
by strings which appear during the QCD phase 
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transition in axion models with iV = 1 assuming 
the axion field is not homogenized by inflation. 
In the simulations, the walls decay immediately, 

i.e. in a time scale of order the light travel time. 
The simulations provide an estimate of the aver¬ 
age energy of the axions emitted in the decay of 
the walls: (wo) — ^rua when '/Xvalma — 20. 

Because of restrictions on the available lattice 
sizes, the simulations are for values of Vajma of 
order 100, whereas in axion models of interest 
Valma is of order 10^®. To address this problem, 
we have investigated the dependence of {u!a)/rna 
upon V^Va/ rua and found that it increases ap¬ 
proximately as the logarithm of this quantity. 
This is because the decay process occurs in two 
steps: 1) wall energy converts into moving string 
energy because the wall accelerates the string, 
and the energy spectrum hardens accordingly; 2) 
the strings annihilate into axions without quali¬ 
tative change in the energy spectrum. If this be¬ 
haviour persists all the way to VXvalma ~ 10^®, 
then {uJa)lrna 60 for axion models of interest. 

This has interesting consequences for the axion 
cosmology and the cavity detector experiments of 
dark matter axion as discussed in Ref. [|| . 
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